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Abstract We processed the data about radial veloci- 
ties and HI linewidths for 1678 flat edge-on spirals from 
the Revised Flat Galaxy Catalogue. We obtained the 
parameters of the multipole components of large-scale 
velocity field of collective non- Hubble galaxy motion as 
well as the parameters of the generalized Tully-Fisher 
relationship in the "HI line width - linear diameter" 
version. All the calculations were performed indepen- 
dently in the framework of three models, where the mul- 
tipole decomposition of the galaxy velocity field was 
limited to a dipole, quadrupole and octopole terms re- 
spectively. We showed that both the quadrupole and 
the octopole components are statistically significant. 

On the basis of the compiled list of peculiar veloc- 
ities of 1623 galaxies we obtained the estimations of 
cosmological parameters rim and cts- This estimation 
is obtained in both graphical form and as a constraint 
of the value 5*8 = {nm/0.3)°-^^as = 0.91 ± 0.05. 

Keyv^rords galaxies; large-scale structure; collective 
motion of galaxies; cosmological parameters 



1 Introduction 

The distribution of matter, including dark matter, in 
the Universe is inhomogeneous on the scales of about 
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100 Mpc. This is manifested, for example, in the ex- 
istence of superclusters and voids. A galaxy, besides 
the cosmological expansion, is also attracted to the re- 
gions with greater density. As a result, the galaxies are 
involved in a non-Hubble large-scale collective motion 
on the background of Hubble expansion. In a more 
sophisticated way this can be considered as a develop- 
ment of initial density and velocity fluctuations in the 
early Universe due to the gravitational instability. 

Investigation of such a motion is important at least 
for two reasons. First of all, it allows to plot the dis- 
tribution of matter, including dark matter, in the sur- 
rounding region of the Universe, which is the main goal 
of cosmography, and to compare this distribution with 
the distribution of luminous matter. The second rea- 
son is that the parameters of this motion are linked 
with certain cosmological parameters, so we can obtain 
new independent estimations of these parameters. Of 
course, the accuracy of such estimation will not be very 
high, but the agreement of different estimations of cos- 
mological parameters can support the correctness of the 
standard c osmological mo del. 

In 1989 iKarachentsevI proposed to use flat edge-on 
spiral galaxies as a tool for studying their large-scale 
collective motion. They are good in this role for the 
following reasons: 

1. The linear diameter is strongly correlated with the 
HI lincwidth for thin bulgeless galaxies. This allows 
to determine distances without photometric data. 

2. Flat galaxies can be easily identifled by their axes 
ratio. 

3. Flat galaxies have a nearly 100 per cent HI detection 
rate. 

4. Galaxies in clusters are usually not flat due to inter- 
action with neighbours. This means that flat galax- 
ies avoid clustering and do not interact with the in- 
tergalactic gas in clusters. 
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5. Peculiar velocities of isolated flat galaxies are not 
perturbed by neighbours. 

To i mplement this idea the F lat Galaxy Catalogue 
(FGC. iKarachentsev et al.llf993l ) was created. It con- 
tained data about N = 4455 galaxies, which satisfied 
the conditions at/bi, > 7 and ai, > 0.6'. Here a?, and bb 
are the major and minor axial diameters directly mea- 
sured on POSS-I and ESO/SERC plates. In accordance 
with the original photographic material, the Catalogue 
consists of two parts: FGC (N — 2573) and its southern 
extension, FGCE {N = 1882). The first part is based 
on the POSS-I and covers the sky region with declina- 
tions between — 20deg and +90deg. The second one is 
based on the ESO/SERC and covers the rest of the sky 
area up to (5 = — 90deg. 

After thorough studies of the catalogue's properties, 
both these parts were joined. The angular diameters 
from the FGCE were converted to the POSS-I sys- 
tem of the FGC, which appeared to be close to the 
a25 system. This system, where galaxy size is taken at 
B = 25 maq/O" isophota l level , was used, in particular, 
bv lde Vaucouleurs et al.l (|l99ll ). Some FGCE galaxies, 
which did not satisfy the condition a > 0.6' after con- 
version, were removed from the resulting Revised Fla t 
Galaxy Catalogue (RFGC, IKarachentsev et aL 1999). 
It contained data about N = 4236 galaxies includ- 
ing the information on the following parameters: Right 
Ascension and Declination for the epochs J2000.0 and 
B1950.0, galactic longitude and latitude, major and mi- 
nor blue and red diameters in arcminutes in the POSS- 
I diameter system, morphological type of the spiral 
galaxies according to the Hubble classification, index 
of the mean surface brightness (I - high, IV - very low) 
and some other parameters, which are not used in this 
article. More detailed description of the catalog ue can 
be found in the paper (jKarachentsev et al.lll999l ). 

The original goal of this catalogue was to estimate 
the distance to galaxies according to the TuUy-Fisher 
relationship in the "HI line width - linear diameter" 
version without using their redshifts. The difference 
between the velocity V derived from the redshift and 
the Hubble velocity R = Hr corresponding to the dis- 
tance r estimated by Tully-Fisher relationship is called 
a peculiar velocity Vpec = cz — Hr. We can use such 
a simple form of the Hubble's law because our sample 
has z < 0.1. 

There are some important things to take into ac- 
count about peculiar velocities. The redshift includes 
not only the velocity of the galaxy, but also the veloci- 
ties of our Galaxy, Solar System and the Earth. Thus, 
to eliminate these factors, all velocities were reduced 
to the frame, where CMB is isotropic. Naturally, the 
redshift gives us only the radial component of the veloc- 
ity and the tangential components cannot be measured. 



Additionally, Tully-Fisher relationship is statistical and 
thus has a certain error. Thus, we can only estimate 
the peculiar velocity for each galaxy, sometimes with a 
significant error. 

However, we believe that there is a large-scale ve- 
locity field. We consider the individual galaxies as test 
particles in the velocity field of large-scale collective mo- 
tion. Thus, having data about the peculiar velocities of 
a large number of galaxies, we can restore the radial 
component of the velocity field. Using some additional 
assumptions, like the potentiality of the flux, it is pos- 
sible to restore the 3D velocity field. For this reason we 
need ample samples of peculiar velocities, preferably 
uniformly covering the celestial sphere. 

These conditions are satisfied by the RFGC cata- 
logue, which covers the whole celestial sphere in both 
hemispheres with the natural exclusion of the Milky 
Way region. However, not all of the RFGC galaxies 
have data required to estimate the distance to them 
using the Tully-Fisher relationship. Nevertheless, the 
sample of galaxies having such data is also quite uni- 
form, as shown on Fig. [T] 

Naturally, RFGC is not the only sample that 
can be used for th is purpose . Possible alterna- 



tives include SBF (ITonrv et a 
ddaCosta et al.' '2000al), SFI ( 
SFI-h+ iSprjiigob ct^ 



1997 



200C), ENEAR 



Havnes e t al.' '1999. 



20071 ). E FAR (Wegner et all 




1996: IColless et al.ll200lh. SMAC (iHudson et al.lll99fl 



|200j) ■ 2MFGC (jMitronova eraP \2Q0^- However 
this article we use only the RFGC catalogue. 

To apply the Tully-Fisher relationship the data from 
the catalogue is not sufficient; we also need to know 
the HI linewidth W (in this article we take it at 50 per 
cent level), or the gas rotation velocity Vrot obtained 
from optical observations. Additionally, we need red- 
shift data. The number of galaxies with such data in- 
creases constantly. 

Since FGC and RFGC were assembled, a num- 
ber of articles were published dealing with collective 
motions of RFGC galaxies. Very preliminary results 



were published by IKarachentsev et al.l (|l995l ). The 



parameters of the bulk m otion were calculated by 
iKarachentsev et al. J 2000b ). In the paper 
(|Parnovskv et al.l 2001 ) not only new data were added, 
but also the models featuring the quadrupole and 
octopole components of the velocity field were pro- 
posed. Also, the generalized Tully-Fisher relationship 
for RFGC galaxies was finalized. It included data not 
only about HI linewidth and angular diameter in red 
and blue bands, but also about the morphological type 
of the galaxy and its surface brightness index. By that 
time the authors had information about radial veloci- 
ties and HI 21 cm line widths or Vrot for 1327 RFGC 
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Fig. 1 Distribution of galaxies over the celestial sphere in galactic coordinates 



galaxies from the difFerent sources listed in the paper 



( Karachentsev et al. 2000bl ). From this number, 1271 



galaxies were included in the sample, and the rest of 
them were considered to be outliers. As a result, the 
first list of peculiar velocities o f RFG C galaxies was pre- 
pared bv lKarachentsev et alj (|2000al ). Four years later, 
the number of gal axies with available data increased 
and reached 1561 |Parnovskv fc Tugavl liool : 1493 of 
them entered the sample. A new version of the list of 



peculi ar velocities was prepared bv lParnovskv &: Tugav 
This list was the basis for solving the two above- 
mentioned problems, namely mapping the matter den- 
sity and estimation of cosmological parameters. The 
distribution of matter densi ty was obtained in the paper 
( Sharov fc Parnovskvlbood ) up to 75h ^ Mpc in the su- 
pergalactical plane and 8 more planes. In the same ar- 
ticle, the excess masses of attractors in this region were 
estimated, and the estimation of /3 parameter was ob- 
tained. The estimations of the cosmol ogical parameters 



flm. a nd (78 were obtained in the paper (jParnovskv et al 



20061) . The model was expanded taking into account 
the effects of general theory of relativity in the paper 
jParnovskv fc Gavdamaka,2004i ). 

An important result , obtained in the paper 
J Parnovskv et a l. '2001") , and confirmed by 
Parnovskv fc Tugav, (2004) . was the statistical signifi- 
cance of quadrupole and octopole components (both at 
more than 99 per cent confidence level). This result 
relied on the F-test, which assumes normal distribu- 



tion of errors (see Hudson 1964h . However, even the 
normally distributed errors of angular diameters and 
HI linewidths and deviations from TuUy-Fisher rela- 
tionship lead to non-Gaussian distribution of peculiar 
velocities. Thus, the statistical significance of these 
ter ms required additional consideration. In the arti- 
cle (jParnovskv fc Parnow ski 2008) it was shown us- 
ing Monte-Carlo simulations that the statistical signif- 
icance of these components appeared to be less than 
perceived from the F-test, but still large enough to be 
considered: the quadrupole was statistically significant 
at 96 per cent confidence level and the octopole - at 
90.. 93 per cent confidence level. 

In the four years that passed since the last list of 
peculiar velocities of RFGC galaxies was compiled, not 
only new data have appeared, but also some old data 
were remeasured. This led to a necessity of repeating 
all the steps required to obtain the new list. This means 
that not only new data have to be added, but the pro- 
cess of sample selection must be repeated from the very 
beginning. As it will be shown further, some results ap- 
peared to be notably different from the ones obtained 
with the previous sample. 

In this paper we present the parameters of the large- 
scale collective velocity field in the framework of the 
multipole model as well as the estimations of some cos- 
mological parameters. The distribution of matter den- 
sity will be discussed in a separate article. 
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2 Multipole models of collective motion of 
galaxies 



In the paper (jParnovskv et al.l 120011 ) the velocity field 
components in the CMB reference frame were expanded 
in terms of galaxy's radial vector r. The first three 
terms are 



Vi = Di + Aikrk + Olf^irkn. 



(1) 



In our notation we use the Einstein rule: summation by 
all the repeating indices. Hence, for the radial velocity 
VsK in the CMB reference frame we get 



(2) 



where n = f/r and r — \r\. Let us decompose the 
tensor A in two parts: Aik — HSik + Q^j.. Here H is 
a trace, corresponding to the Hubble constant, Sik is 
the Kronecker delta, and Q' is a traceless tensor. Next 
we switch from distance r to the corresponding Hubble 
velocity R — Hr. 



VsK = R + DiUi + RQikUink + R OikiniUkni, 



(3) 



where Q = H-^Q' and O = R-^O'. 

This decomposition was used to obtain some mod- 
els of dependence of galaxy's radial velocity V. In the 
simplest D-model (Hubble law + dipole) we have 



V3K^R + V^'P + SV, V^'P = A 



(4) 



where D is a velocity of homogeneous bulk motion, SV 
is a random deviation and n is a unit vector towards ga- 
laxy. After the addition of quadrupole terms we obtain 
a DQ-model 

V3K=R + V^'P + + 6V, V^''"^ = RQ^kU^Uk (5) 

with symmetrical traceless tensor Q describing quadru- 
pole components of the velocity field. It characterises 
the relative deviation of an effective "Hubble constant" 
in a given direction from the mean value. More de- 
tailed discussion will be given further. The DQO-model 
includes the octopole component of velocity field de- 
scribed by a symmetrical tensor O of rank 3: 



V^ii = R + V^'P + V"^" 
y-* = R^OMn,nkni. 



6V, 



(6) 



In some cases it makes sense to use another way of 
describing the octopole component in the DQO-model. 
For this purpose the tensor O can be considered a sum 
of a traceless tensor O, and a trace, characterized by 
the vector P 



ijk 



o 



ijk 



ijkOjk- 



(7) 



Here the indices in parantheses are symmetrized. Thus 
an alternative form of the DQO-model is 



V3K = R+ (A + R'^Pt)n^ + RQ.kn^nk 
+ R^dikiniUkni -f 5V. 



(8) 



We will use Cartesian components of the vector n in 
the galactic coordinates: 



ni = Tiz = sin h, n2 = rix = cos I cos 5, 
= fT-y = sin I cos b. 



(9) 



These three models were used for processing data on 
RFGC galaxies. To estimate the distances to galaxies 
a generalized TuUy-Fisher relationship was used in the 
'linear diameter - HI line width' version. It has a form 



( Parnovskv et al.ll200lh 



R={Ci+ C2B 



C3BT 



(10) 



where is a corrected HI line width in kms~^ mea- 
sured at 50 per cent of the maximum, a is a cor- 
rected major galaxies' angular diameter in arcminutes 
on red POSS and ESO/SERC reproductions, [/ is a 
ratio of major galaxies' angular diameters on red and 
blue reproductions, T is a morphological type indicator 
{T — It — 5.35, where It is a Hubble type; It = 5 
corresponds to type Sc), and i? is a surface bright- 
ness indicator {B — Isb — 2, where Isb is a surface 
brightness index from RFGC; brightness decreases from 
I to IV) . Details of corrections one can find in the pa - 



pers ( Parnovskv et al.ll2001 : Parnovskv fc Tuga'vl2004 ). 
The reasons for choosing this form of Tully-Fisher re- 
lationship for RFGC galaxies are given ibid. We only 
note that the statistical significance of each term in eq. 
(fTO|) is greater than 99 per cent. 

Note that eq. (fTO|) differs from the classical Tully- 
Fisher relationship, which has the form 



R^C 



(11) 



The corresponding ana lysis was performed in the pa- 
per (jParnovskv fc Ga vdamaka 2004) for the relativis- 
tic model of motion. It was shown that the power a 
for the classical form of Tully-Fisher relationship (fTT|) 
somewhat differs from 1. However, after including the 
corrections for the surface brightness and the morpho- 
logical type of the galaxy, the value of a appeared to 
be close enough to 1 to use the form given by (fTU|) . If 
we repeat this analysis for the new data, the minimal 
variance is reached at a ~ 1.16. At the same time the 
standard deviation a for the form (fTTj) is 5.. 7 per cent 
more than for the form (jlOp depending on the depth 
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of the subsample and model used. If we combine the 
forms (Uni) and lU) into 



R^{C\ + C2B 



C3BT + CaU) 



(12) 



we obtain the minimal variance at a ^ 1.05, which is 
very close to 1. This correction does not reduce a due 
to the decreased number of degrees of freedom. This 
makes the introduction of a senseless. For this reason 
we will use only the form (fTO|) . This is additionally 
justified by the fact that all the coefficients enter this 
equation linearly and thus we can use the linear regres- 
sion analysis. 

The D-model has 9 parameters (6 coefficients C and 
3 components of vector D), DQ-model has 14 parame- 
ters (5 components of tensor Q are added), and DQO- 
model is described by 24 coefficients. The values and 
errors of the coefficients were calculated by the least 
square method for different subsamples with distances 
limitation to make the sarnple rn ore homogeneous in 
the depth (jParnovskv et al.l 120011 ). We used the same 
weights for all datapoints. Since the quadrupole and 
octopole terms explicitly depend on i?, an iteration pro- 
cedure was used. Note that the coefficients of the gener- 
alised TuUy-Fisher relationship and components of the 
velocity model were fitted simultaneously. The itera- 
tions converge rather quickly. 



3 Description of data 

Directly from the RFGC catalogue we take the an- 
gular diameters on blue and red reproductions, mor- 
phological type of a galaxy and its surface brightness 
index. The angular diameters are corrected for in- 
trinsic absorption. From other sources we take the 
radial velocity in the CMB frame and the HI line 
width or the gas rotation velocity Vrot obtained from 
optical observations. The line widths are corrected 
for cosmological expansion and for turbulence. The 
data on HI line widths fo r RF GC galaxies in the 



HyperLeda (jPaturel et al.l l2003al ) catalogue are based 
on 149 sources, which have different reliability. The 
main sources used in previou s versions of our sam- 



ples are listed in the papers (iParnovskv et al.l 12001 



Karachentsev et al 



2000ai b 



Parnovskv fc T ugav 200. 
Parnovskv fc Tugav . 20051 ). Since that time new large 
volumes of data on radial velocities and HI linewidths 
of galaxies were published (e.g., ?). For this reason we 
rechecked the original data on HI line widths of galax- 
ies in previous samples and changed some of them. For 
each galaxy we used only the data of original measure- 
ments and did not average the data of different papers. 



In the cases when data from reliable sources differed 
significantly, we used less reliable sources to choose be- 
tween them. We also checked the deviation of the dis- 
tance to galaxy from the TuUy-Fisher relationship in 
the D-model. 

After revising data from the previous samples, we 
added new data measured during the last 5 years. 
We took HI linewidths from the original articles 



i Koribalski et al.l2 004VTh eureau et al.ll2005HDovle et al.| 
2005t ISpringob et al.ii2005f ). In the case when W was 
unavailable we converted from the HyperLeda cat- 



alogue (jPaturel et al.ll2003bl ) to W using the relation- 
ship W = 2.016 * Vrot + 20.6, obtained by processing 
data for galaxies, which have both records. Note that 
there are rare cases of significant deviations from this 
relationship. 

In the process of data preparation we constantly 
checked the deviations of the radial velocities from the 
velocities calculated in the D-model. Some data were 
obvious outliers. Such galaxies were rejected from the 
sample. As a result, we had all the necessary data for 
1678 galaxies, including 55 outliers. Thus the sample 
contained 1623 galaxies, including 204 new and 1419 
old ones. It is noteworthy that some of the old data 
changed as well. The distribution of these galaxies over 
the celestial sphere is shown on Fig. [TJ The mean ra- 
dial velocity of old galaxies is 6000 km s^^, and of new 
ones it is 5600 km s~^. The mode in both cases is about 
4800 /cm s~^. The sample has 90% completeness up to 
about 6000 fcms"^. 

The farthest galaxies in the sample have radial veloc- 
ities above 20000 fcrns^^. However, at such scales the 
sample is inhomogeneous and incomplete. To deal with 
more homogeneous samples, we consider subsamples 
with distance limitation. They include all the galaxies 
satisfying the condition R < Rmax in the D-model and 
are used for all the three models. For this purpose we 
used the following iterative procedure. First, we calcu- 
late the distances to the galaxies using seed coefficients 
from the whole sample or another subsample with close 
Rmax- Then we select the galaxies with R < Rmax and 
recalculate the coefficients using only their data, once 
again calculate the distances and so on. In most cases 
these iterations converge, but in some cases, more of- 
ten for small Rmax, we get a limit cycle. In this case 
we select only the galaxies whose distances are below 
Rmax for all sets of coefficients in this limit cycle. 

We used different subsamples with Rmax rang- 
ing from 3000 to 12000 fcm s~^. In this paper we 
mostly present information for the subsamples with 
Rmax equal to 8000 fcms^^ (r = 80 ^ Mpc) and 
10000 fcms-i (r ^ lOOh^^ Mpc). 
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Fig. 2 Dependence of the dipole component of the velocity 
on the sample depth for 3 models. Triangles correspond to 
the D-model, squares - to the DQ-model, and circles - to 
the DQO-model 
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Fig. 3 Bulk motion apices in galactic coordinates for 
Rmax = 10000 fcms~^. Crosses mark the apices of the bulk 
motion in the D-model surrounded by la, 2a and Scr con- 
fidence areas. The square marks the apex position of the 
galaxies from the previous sample with corrected data. The 
diamond marks the apex position in the DQO-model. Solid 
boundaries corr espond to new result s , and dashed ones - 
to the results of IParnovskv fc Tugavl (|2004|). Numbers de- 
note the results of other authors: 1 - IILvnden-Bell et al 



1988h . 2 - jHudson et all I" 
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1994f ) 4 - (IParnovskv et all 



Ida Costa et al. 



Dale et a 



IWatkinsetal .1120091 ) 



995h 



3 

200 ll ). 5 



Lauer fc Postman 



llDekel et al.l 



2000b|), 7 - jHudson et al I 
1999l l. 9 - (|Kudrva et all l2003h . 



199 



200. 
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4 Multipole structure of collective motion 

The parameters of a generalized Tully-Fisher relation- 
ship for all three models and two Rmax values are given 
in Table [T] together with their statistical significance 
levels according to F-test. The F values should be com- 
pared to 3.8, 6.6, 7.9, 10.8 and 12.1, which correspond 
to 95, 99, 99.5, 99.9 and 99.95 per cent confidence levels 
respectively. Most of them are statistically significant 
at 99.95 per cent confidence level with the exception 
of Cs for Rmax = 1 0000 fcms~^. Comparing the coef- 



ficients to that of IParnovskv fc Tugavl (|2004l ) one can 
see that Ci, C4, C5 and Cq remain the same within the 
margin of error, and C2 and especially C3 are reduced, 
which means lesser influence of the morphological type 
of galaxies. This also leads to some decrease of statis- 
tical significance, which is still more than 95 per cent. 

The inclusion of the octopole component for sub- 
samples with Rmax about 8000 fcms"^ leads to a 
drastic d ecrease of the d i pole c omponent up to Icr 
level, 



pie, which lead to dipole and octopole components be- 
ing non-orthogonal. 

This effect is important because D and DQO models 
yield the values of the bulk motion velocity differing by 
a factor of 2. Thus the dependence of the bulk motion 
velocity on the sample depth will be totally different 
for D and DQO models as one can see from Fig. [51 
One should take this into account when analysing the 
convergence scale. 



In com parison with the papers (|Parnovskv et al 



20011 ) and (jParnovskv k, Tugavl 120041 ) . the apex direc- 
tion of the dipole component has changed. In Table 
[5] we present its parameters for 2 subsamples together 
with the standard deviation a for D-model. Note that 
a includes both the intrinsic scatter in the generalised 
Tully-Fisher relationship and the stochastic component 
of the velocities, resulting in high values. 

The comparison o f the results from Table H] with 

those obtained earlier I Parnovskv et al.ll2001 : Parnovskv fc Tugav^ 



as 



Parnovskv et al 



20041) shows that the norm of the velocity hasn't 



( 20011 ) have shown. This 
effec t appeared to be even rn ore evident in the pa- 
per (jParnovskv fc Tugavl |2004[ ). Such a strong dipole- 
octopole coupling can be explained by the same sym- 
metry of these components with respect to the inver- 
sion of space. For the new sample this effect appeared 
to be weaker as the dipole component decreases only 
up to 2(7 level (see Fig. [5]). Naturally this effect is 
caused by incompleteness and asymmetry of the sam- direction has changed at slightly more than 2cf. Besides, 



changed significantly, but the apex direction has. Let us 
discuss this question in detail. On Fig. |3]we present the 
MoUweide projection upon the celestial sphere of the el- 
lipsoids corresponding to Icr, 2a and Scr confidence lev- 
els for the apex direction in the framework of D-model 
for Rmax = 10000 km s~^. Solid boundaries correspond 
to new results, and dashed ones - to the results of 



Parnovskv fc Tugavl (|2004l ) . One can see that the apex 
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Table 1 Parameters of the generalized TuUy-Fisher relationship (|10p 



% 


D-model 


DQ-model 


DQO-model 


c. 








Cl 


Fi 








Rrnax = 8000 km/s 


(1240 


galaxies) 








1 


15.8 ± 


1.6 


103 


17.2 ± 1.6 






122 


17 A ± 


1.6 


126 


2 


1.84 ± 


0.20 


80.5 


1.96 ± 0.21 






90.0 


1.95 ± 


0.21 


88.4 


3 


-0.46 ± 


0.12 


13.4 


-0.49 ± 0.13 






15.2 


-0.45 ± 


0.13 


12.9 


4 


8.29 ± 


1.33 


38.6 


7.13 ± 1.34 






28.6 


7.24 ± 


1.33 


29.6 


5 


(-5.54 ± 


1.53) X 10-3 


13.1 


(-6.63 ± 1.53) 


X 10- 


3 


18.7 


(-7.85 ± 


1.57) X 10-3 


25.0 


6 


-825 ± 


86 


92.4 


-867 ± 85 






103 


-883 ± 


85 


107 








Rmax = 10000 km/s 


(1459 


galaxies) 








1 


16.3 ± 


1.4 


143 


17.3 ± 1.4 






157 


17.2 ± 


1.4 


156 


2 


1.38 ± 


0.18 


59.3 


1.44 ± 0.18 






62.1 


1.55 ± 


0.18 


70.6 


3 


-0.24 ± 


0.11 


4.55 


-0.26 ± 0.11 






5.32 


-0.27 ± 


0.11 


5.54 


4 


7.62 ± 


1.17 


42.4 


6.80 ± 1.18 






33.0 


6.99 ± 


1.18 


35.3 


5 


(-5.96 ± 


1.05) X 10-3 


32.0 


(-6.41 ± 1.06) 


X 10- 


3 


36.4 


(-6.71 ± 


1.13) X 10-3 


34.9 


6 


-785 ± 


87 


81.8 


-817 ± 87 






88.2 


-845 ± 


87 


94.9 



Table 2 Parameters of the dipole component 



Model 


Rmax 




Dx 


Dy 


D 


I, deg 


b, deg 


fj, km/s 


km/s 


D 


8000 


147 ± 47 


176 ± 58 


-136 ± 56 


267 ± 65 


322 


33 


1045 


DQ 


8000 


81 ± 50 


174 ± 60 


-134 ± 58 


234 ± 52 


323 


20 


1033 


DQO 


8000 


-59 ± 79 


122 ± 104 


-103 ± 100 


170 ± 95 


320 


-20 


1020 


D 


10000 


139 ± 49 


229 ± 59 


-106 ± 59 


288 ± 64 


335 


29 


1159 


DQ 


10000 


70 ± 53 


200 ± 62 


-116 ± 62 


242 ± 56 


330 


17 


1151 


DQO 


10000 


61 ± 77 


85 ± 98 


-182 ± 97 


210 ± 86 


295 


17 


1138 



on the Fig [HI we show th e resu lts of other authors : 1 



( Lvnden-Bell et al 



Lauer fc: Postma: 



Dekel et al 



i 



1988), 2- 
19941) ■ 4 
1999h. 6 - 



(Hudson et al 



iudson et al.lll995l). 3 



(Ida Costa et all l2000bD 



2004h . 8 - dOale et al 



Parnovskv et al l l200 ih . 



1998 



j_ — . — 9 - 

( Kudrva et al.l 120031 ). 10 - (IWatkins et all 12003) One 
that the result of IParnovskv fc TukotI \2004) 



can see 
matches other results better. 



The new result became 



closer to the result of lLauer fc PostmanI (|l994f ) , but our 
norm is 2 times less. What caused such a deviation? 
We found the apex for 'old' galaxies used in previous 
samples but with corrected data for the HI lincwidth. 
At this stage the apex has notably deviated from its 
original position. This position is shown on Fig. [3] by a 
square. After adding new data this deviation increased. 

Let us discuss how much this deviation affects the 
global picture of the large-scale motion. Some au- 
thors believe that the norm and the apex of the di- 
pole component are the most prominent characteristics 
of large-scale motion. This was true for very old sam- 
ples with small depth. Our sample, however, has so 
large depth that it includes most nearest superclasters 
like the Great Attractor and Perseus-Pisces. For this 



reason, the first few components of the multipole de- 
composition ^ can not adequately describe all the ra- 
dial velocity field of the large-scale motion. Moreso, 
this is hard to do with a single component This 
phenomenon is due to the fact that a multipole decom- 
position is capable of describing very-large-scale motion 
with characteristic scale much larger than the scale of 
the sample, but can not describe the not-so-large-scale 
motions with characteristic scale less than the scale of 
the sample caused by attraction to these superclusters 
for any reasonable number of components. 

As a result, the variation of the subsample depth 
Rmax can lead to significant deviation of the vector D, 
especially when the supsample boundary crosses a large 
attractor. At smaller Rmax values we observe only the 
increase of radial velocity in the direction of the attrac- 
tor. At larger Rmax we see also the backfall. In the 
framework of DQO-model these effects can be taken 
into account, but it is too much for the simple D-model 
to handle. On Fig. H we show the variation of the apex 
position in dependence from the sample depth. One 
can see that at small Rmax the apex crosses the galac- 
tic plane and approaches the results of most authors 
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Table 3 Coefficients of the quadrupole component 



i 


Rmax = 8000 km/s 


Rrnax = 10000 km/s 


DQ-model 


DQO-model 


DQ-model 


DQO-model 


q^, % 


F^ 


q^, % 




q^, % 




ft, % 


F^ 


1 


8.1 ± 1.5 


28.8 


7.7 ± 1.6 


22.6 


5.7 ± 1.4 


16.6 


7.6 ± 1.6 


22.6 


2 


-2.1 ± 1.6 


1.6 


-1.9 ± 1.7 


1.2 


-1.9 ± 1.5 


1.7 


-2.5 ± 1.6 


2.3 


3 


-0.8 ± 2.2 


0.1 


-1.1 ± 2.3 


0.2 


2.1 ± 1.9 


1.2 


-1.2 ± 2.1 


0.3 


4 


0.9 ± 2.3 


0.1 


1.2 ± 2.5 


0.2 


4.1 ± 2.1 


4.0 


2.6 ± 2.5 


1.1 


5 


1.7 ± 2.7 


0.4 


0.1 ± 2.8 


0.0 


1.7 ± 2.5 


0.5 


-0.2 ± 2.7 


0.0 


Total 




6.44 




5.10 




4.63 




4.85 



Table 4 Eigenvalues and eigenvectors of the quadrupole tensor 



Rmax ; 


N 


Maximum 


Minimum 


Third axis 


km/s 


Qi, % 


I 


b 


Q2, % 


I 


b 


Qs, % 


I 


b 



DQ-model 



3000 


292 


19 


2 


± 


3 


7 


151 


-7 


-19.8 ± 4.2 


80 


68 





6 


± 


5 


6 


58 


-20 


6000 


916 


5 


6 


± 


1 


6 


186 


74 


-8.5 ± 1.8 


59 


10 


2 


9 


± 


2 


4 


327 


13 


7000 


1068 


7 


3 


± 


1 


5 


198 


81 


-5.9 ± 1.7 


77 


4 


-1 


3 


± 


2 


3 


347 


7 


8000 


1240 


8 


1 


± 


1 


4 


142 


87 


-6.2 ± 1.5 


102 


-2 


-1 


9 


± 


2 


1 


12 


2 


8500 


1322 


6 


8 


± 


1 


3 


68 


85 


-5.2 ± 1.5 


106 


-4 


-1 


6 


± 


2 





16 


-3 


9000 


1379 


7 





± 


1 


2 


53 


79 


-5.6 ± 1.4 


109 


-6 


-1 


4 


± 


1 


8 


18 


-9 


9500 


1417 


6 


9 


± 


1 


2 


41 


78 


-5.4 ± 1.4 


109 


-4 


-1 


5 


± 


1 


8 


18 


-11 


10000 


1459 


6 


3 


± 


1 


2 


55 


75 


-4.4 ± 1.3 
DQO-model 


105 


-10 


-1 


9 


± 


1 


8 


13 


-11 


3000 


292 


22 


3 


± 


4 


3 


153 


-9 


-19.9 ± 4.3 


80 


63 


-2 


4 


± 


6 


1 


59 


-26 


6000 


916 


5 


1 


± 


1 


9 


156 


71 


-7.8 ± 1.9 


56 


3 


2 


8 


± 


2 


7 


325 


19 


7000 


1068 


6 


6 


± 


1 


8 


124 


81 


-5.4 ± 1.7 


75 


-6 


-1 


2 


± 


2 


5 


346 


7 


8000 


1240 


7 


8 


± 


1 


6 


141 


86 


-5.9 ± 1.6 


90 


-3 


-1 


9 


± 


2 


3 


360 


3 


8500 


1322 


7 


5 


± 


1 


6 


133 


84 


-4.9 ± 1.5 


92 


-5 


-2 


6 


± 


2 


2 


2 


4 


9000 


1379 


8 


2 


± 


1 


6 


101 


84 


-5.1 ± 1.5 


90 


-6 


-3 


1 


± 


2 


2 


359 


1 


9500 


1417 


7 


8 


± 


1 


6 


93 


86 


-4.8 ± 1.5 


99 


-4 


-3 





± 


2 


2 


9 


-1 


10000 


1459 


7 


8 


± 


1 


6 


120 


83 


-5.3 ± 1.5 


89 


-6 


-2 


5 


± 


2 


2 


359 


3 



depicted on Fig. [3] An additional argument in favour 
of the apex deviation being not much of a problem is the 
following consideration. If in the framework of D-model 
with the same Rmax for each galaxy we introduce statis- 
tical weights corresponding to the variance of the radial 
velocity as a function of distance, the apex will corre- 
spond to the position I = 327 deg; b = 18deg, w hich is 



closer to the result of iParnovskv fc Tugavl (|2004f ). It is 
obvious that in such an approach the input from the 
nearest galaxies will dominate over the more distant 
galaxies. 

Another possible explanation is that the velocity 
field is not well described by a pure dipole, and the 
dipole one recovers is unstable with respect to the limit- 
ing distance of the sample used. Furthermore, although 
the dipole is (in principle) a well-defined property of the 
field, and its mean convergence with distance is pre- 
dictable within a given cosmological model, the prob- 



lem here (as noted previously) is that the multipole ba- 
sis is not orthogonal for this sample, so that in practice 
the dipole is not well-defined (in addit ion to varyin; 



stoch astically with sample volume) . As (jWatkins et al.l 
20091) have shown, this sample bias in the measurement 



of the dipole is substantial, and has to be taken into ac- 
count when comparing results from different samples. 

As we mentioned before, the usage of the more 
complex DQO-model allows to describe such phenom- 
ena as the backfall. Beside s, as it was shown by 
Parnovskv fc Parnowski ( 20081 ). the mean deviation of 



the calculated apex from its true value, caused by mea- 
surement errors and the deviations from the TuUy- 
Fisher relationship, is less in the DQO-model in com- 
parison to the D-model. On Fig. [3] the relevant apex 
position for the DQO-model is denoted by a diamond. 
We can see that it fairly corresponds to the results of 
other authors. 
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Table 5 Coefficients of the octopole component 



-90° 



Fig. 4 Dependence of bullc motion apex on tlie sample 
depth 



As it was mentioned above, we used the generalised 
Tully-Fisher relationship in the form 1101 However, we 
also calculated the bulk motion parameters in the D- 
model using the classical form[TlJ We obtained that the 
bulk motion velocity is 274 fcms~^ and the apex coor- 
dinates are I = 317, 5 = 36 for Rmax — 10000 fcms"^. 
For Rmax — 8000 fcms^^ the corresponding values are 
254 fcms^^, I — 306, b = 41. As one can see from Table 
[Hand Fig. [3] these values fall in Icr confidence areas. 

Let us consider the quadrupole component. In Table 
[3] we present the corresponding coefficients. The trace- 
less quadrupole tensor is described by 5 coefficients: 



V;?'- = RQ,jn,nj = R[qi{n\ - nl) + 92 ("1 
- nl) + q3nin2 + g4"-i'^3 + g5n-2'^3 



(13) 



In Table [3] we also present the errors of these coeffi- 
cients and their significances according to F-test. The 
quadrupole coefficients Qj are close to the coefficients 
obtained by Parnovskv &: Tugav ( 2004h . One can see 
that the significance of all components except the first 
one is low. Nevertheless, the significance of the first 
component is high enough for the total significance of 
the whole quadrupole component to be sufficient. The 
F values of the quadrupole component are the follow- 
ing: 4.63 for the DQ-model and Rmax — 10000 fcms"^, 
4.85 for the DQO-model and Rmax = lOOOO km s~'^, 
6.44 for the DQ-model and Rmax = 8000 fcms~^, and 
5.10 for the DQO-model and Rmax = 8000 km s''^. 
This should be compared to 3.35, which corresponds 



1 


Rrnax 


= 8000 km/s 


Rmax 


= 10000 km/s 




10^ 






■ 10^ 


F^ 


1 


7.2 


± 3.9 


3.5 


-4.1 


± 3.0 


1.8 


2 


-2.1 


± 5.2 


0.2 


-0.4 


± 3.7 


0.0 


3 


-3.7 


± 5.4 


0.5 


-1.6 


± 4.6 


0.1 


4 


20.3 


± 7.1 


8.2 


10.6 


± 4.9 


4.7 


5 


-7.8 


± 7.4 


1.1 


4.3 


± 6.1 


0.5 


6 


24.0 


± 7.6 


9.8 


23.0 


± 5.5 


17.5 


7 


-13.0 


± 9.6 


1.8 


-12.5 


± 8.0 


2.4 


8 


10.9 


± 7.8 


2.0 


15.8 


± 6.1 


6.8 


9 


-17.9 


± 9.0 


3.9 


-6.2 


± 6.9 


0.8 


10 


33.8 


± 10.2 


11.1 


25.9 


± 7.6 


11.7 


Total 






4.16 






4.37 



to the confidence level of 99.5 per cent for 5 degrees 
of freedom. Thus the quadrupole component is statis- 
tically significant according to F-test at 99.5 per cent 
confidence level both in DQ and DQO-models. 

What is the physical sense of the qu adrupole compo- 



nent? As one can see from the paper (jParnovskv et al 



20011 ). it can be naturally combined with the Hubble 
constant. As a result, we obtain the effective 'Hubble 
constant' depending on direction 



H{l,b) = H{l + Q,kn,nk). 



(14) 



Naturally, this effective 'Hubble constant' is caused by 
large-scale collective motion on the sample scale. To 
estimate the value of its anisotropy we found the eigen- 
values and eigenvectors of tensor Q. The three eigen- 
vectors are orthogonal and the sum of three eigenvalues 
is equal to zero because Q is a traceless tensor. 

In Table 2] we present the maximum, minimum and 
medium eigenvalues and the corresponding eigenvectors 
for DQ and DQO-models and different Rmax- To ob- 
tain the eigenvalues of this 'Hubble constant' we should 
add 100 per cent to each of the listed values. For exam- 



ple, for the subsample with R„ 



10000 km in 



the framework of DQ-model the maximal 'Hubble con- 
stant' is equal to 106.3 per cent of the actual Hubble 
constant, and the minimal - 95.6 per cent. We can see 
that the anisotropy is weak but statistically significant. 
The anisotropy is large only for very small scales less 
than 60h~^ Mpc. As one can see from the Table HI the 
maximum axis is dominant over the other two. How- 
ever, the direction of the minimum axis is more stable 
for both DQ- and DQO-model. In the DQO-model the 
maximum axis is also stable at Rmax > 6000 fcms~^. 
It lies almost in the supergalactic plane. The minimum 
direction deviates from this plane by approximately 30 
degrees of arc. 
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Fig. 5 Distribution of radial velocity of collective non-Hubble motion at 20h ^ Mpc 



In the DQ-model the situation is different. As one 
can see from the Table |4j the variation of the maximum 
eigenvector is more evident. With the increase of Rmax 
the eigenvector starts deviating from the supergalac- 
tic plane and at Rmax = 10000 A:ms~^ this deviation 
reaches approximately 20 degrees of arc. This effect can 
be explained by the fact that the Local Supercluster in 
this case constitutes the smaller fraction of the sample. 

We also calculated the statistical significance of the 
octopole component in the DQO-modcl. The cor- 
responding F values are equal to 4.37 for Rmax — 
10000 fcms^i and 4.16 for Rmax = 8000 /cm s"^ 
These values should be compared to the value 2.52, 
which corresponds to the confidence level of 99.5 per 
cent for 10 degrees of freedom. Thus, the octo- 
pole component is also statistically significant accord- 
ing to F-test at 99.5 per cent confidence level. The 
same qualitative conclusion wa s achieved in the paper 
( Parnovskv &: Parnowskil 2008 ). We can also calculate 
the statistical significance of P in eq. Q as in the p aper 
( Parnovskv et ahlboOll : IParnovskv fc Tugavll2004l) . Its 
F values are equal to 0.21 for Rmax = 10000 km and 
2.50 for Rmax — 8000 km s^^. This should be compared 
to 3.78, 2.60, 2.14 and 0.79, which correspond respec- 
tively to 99, 95, 90 and 50 per cent confidence levels for 
3 degrees of freedom. Thus, for Rmax — 10000 fcms"^ 
the octopole trace is insignificant at 50 per cent confi- 
dence level, and for Rmax = 8000 km it is significant 
at 90 per cent confidence level, but this is insufficient 
to claim that this value is significant for all subsamples. 



In Table [5] we present the coefficients of the octopole 
component, their errors and significances according to 
F-test. The octopole tensor including the trace P is 
described by 10 coefficients: 



= R^dijkniUjnk = R^{din\ + d2n\ 

+ djn2nl + dgn^nl + dgUzn^ 
+ diQnin2n3). 



(15) 



The coefficients of the octopole coni p onent differ 
from those given by Parnovskv &: Tugav ( 2004J l much 
stronger than the quadrupole ones. Nevertheless, their 
la confidence areas overlap. 

The knowledge of the coefficients of dipole, quadru- 
pole and octopole components allows us to calculate 
the distribution of the radial velocity of the collective 
motion of galaxies. On Fig. [5] we present such a distri- 
bution for the DQO-model with Rmax — 10000 fcms^^ 
and R = 2000 fcms"^, and on Fig. [S] - for Rmax = 
10000 fcms~i and R = 8000 km s'^. Note that the am- 
plitudes of the dipole, quadrupole and octopole com- 
ponents are close to each other at i? 6000 fcms^^. 
Thus, for larger distance like in Fig. [S] the octopole 
component prevails in the velocity field. This figure 
is characterised by a deep minimum of radial velocity 
in the direction opposite to the Great Attractor and 
Shapley concentration. In the direction to these super- 
clusters there is also a minimum caused by a backfall 
to the Great Attractor. For smaller distances R like 
in Fig. [5] this minimum becomes a maximum, but the 
mininmm in the opposite direction is deeper than this 
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Fig. 6 Distribution of radial velocity of collective non-Hubble motion at 80h ^ Mpc 



maximum. It is clear that this picture is a simpHfied 
representation of the velocity field, but still it is much 
more complex than the simple bulk motion given by the 
D-model. 



5 Estimation of some cosmological parameters 
and their combinations 

As a result, we obtain not only the parameters of the 
multipole models of the collective motion, but also an 
estimation of the distance to galaxies using the gen- 
eralised TuUy-Fisher relationship. This allows us to 
compile a list of peculiar velocities for 1623 galaxies. 
It can be used to estimate the cosmological parame- 
ters ilm and fjg ■ For the previous ve rsion of the s ample 
it was done in the paper (fParnovskv et al 



use a method similar to that used by iFeldman et al 



2nnm. We 



12003) but with some changes. We calculate a 2-point 
correlation function for peculiar velo cities, which is ap- 



JTOximated by a formula obtained bv lJuszkiewicz et al 
19991 . The formula has a nonlinear dependency on 
and as , thus these parameters can be estimated by the 
maximum likelihood method using the correlation func- 
tion for RFGC galaxies. This yields an estimation of 
cosmological parameters Q,m and erg depicted on Fig. 
[7] There are 2 versions of the approximation formula: 
a simple one and a complex one. For the complex for- 
mula we obtain the values Qm ~ 0.26 and trg = 0.96. 
This point is surrounded by 1(t and 2a confidence ar- 
eas. They have the banana-like shapes - very long and 



narrow curved strips. Judging from the boundaries of 
these areas, the errors of the cosmological parameters 
are very large. For this reason, the good agreement 
in the value of flm between our estimation and other 
estimations can be just incidential. The simple model 
yields similar bananas. 

Let us compare the obtained confidence areas with 
more accurate estimations. Let us start from the 3-year 
WMAP resul ts. On Fig. [5] we reproduce Fig. 7 from 



the paper ([Spergel et al.l 120071 ) with an overlaid plot 



of la and 2a confidence areas for the complex model. 
The blue colour depicts WMAP data, the yellow colour 
depicts data obtained from weak gravitational lensing, 
the green colour depicts constraints imposed by both 
WMAP and lensing, and the red colour depicts our 
constraints. We can see that our results are similar 
to the weak lensing data much more than to WMAP 
data. In any case, these results are obtained with to- 
tally different methods, and their good correspondence 
is very promising. 

The results of comparison with 5-year WMAP data 
are presented o n Fig. M It is a reprod uction of Fig. 19 
from the paper (jKomatsu et al.ll2008l ) with an overlaid 
plot of 1(7 and 2a areas for the complex model. It con- 
tains constraints set by WMAP, mutual constraints set 
by WMAP, BAO and supernovae, and our constraints. 

The main drawback of our estimation is its graph- 
ical form. In many cases it is preferable to deal with 
numerical constraints. To obtain such constraints we 
will make use of the long and narrow shape of our con- 
fidence areas. If we consider a combination of cosmo- 
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Fig. 7 la and 2a confidence levels of cosmological con- 
stants flm and as for the complex (solid) and simple 
(dashed) models 



logical parameters of the form (f2m/0.3)"cr8 and a in 
the range from 0.3 to 0.7 we will obtain an estimation 
of the combination with better accuracy then either of 
the two cosmological parameters. The corresponding 
estimations are given on Fig. [TOl 

For the previous version of the sample the estima- 
tion of such a c ombination was made in the paper 
(|Parnovskvll2008l ). Ibid there are many estimations of 
other authors for different values of a. Of all the val- 
ues of a the smallest error was reached at a = 0.35. 
On Fig. [Tl] one can see the boundaries of la confi- 
dence areas for the complex and simple models in the 
(Om/0.3)° '^^cr8 against ilm coordinates. Both of them 



fit the constraint {nm/0.3)°-^^as = 0.91 ±0.05. This re- 
sult is almost the same as for the previous version of the 
sample (0.93 ±0.05). The importance of this a value is 
caused by the fact tha t the value = (^ m/0■•i)°'^^aH 
was used in the paper ( Evrard et al. 20081 ). As the au- 
thor notes, the low estimation 5*8 = 0.69 based upon 
WMAP and SDSS galaxies clusterisation leads to a 
number of problems. It leads to a contradiction with re- 
sults of clusters formation modelling. The high estima- 
tion 5*8 ~ 0.8 — 0.9 lifts these problems. The estimation 
of the share of hot gas following from it corresponds 
to modern data obtained from the Sunyaev-Zel'dovich 
effect. Our res ults support the h igh estimation. Note 
that the paper ( Evrard et al.|[200 8) contains many es- 
timations obtained by different authors with different 
methods. Among them there are both high estimations, 
close to ours, and low estimations. 
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Fig. 8 An estimation of cosmological constants Qrn and 
as superimposed on 3-year WMAP results 



6 Conclusion 

We prepared a new increased and largely revised sam- 
ple of RFGC galaxies with data about redshifts and HI 
linewidths. It allowed us to improve the estimations 
of parameters of the radial velocity field of large-scale 
collective motions of galaxies using 3 models of its mul- 
tipole structure. In comparison with the previous ver- 
sions the main features remained intact, but separate 
details, for example the apex of the dipole component of 
the bulk motion in the D-model, significantly changed. 
As before, the statistical significance according to F- 
test of both the quadrupole and the octopole compo- 
nents is well above 99.5 per cent. The exact parameters 
are given in the article for subsamples with maximum 
distances lOOh^^ Mpc and SOh'^ Mpc. 

The obtained velocity of the bulk motion is in 
agreement with the ACDM model expectation of ^ 
200 km s^^. The value of the bulk motion velocity 
at tracted additio n al att ention after the recent paper 
of IWatkins et al.l (|2009l ) who obtained this value as 
407 ± 81 km at the same scale 100 h^^ Mpc as in 
present article. In this connection some authors imme- 
diately started speculating about the challenge to the 
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Fig. 9 An estimation of cosmological constants Qm and 
(T8 superimposed on 5-year WMAP results 



ACDM model and the necessity for in clusion of non- 
gravitational forces (lAvaita et al. 20091 ) . 

For 1623 galaxies we compiled a list of peculiar ve- 
locities which will be published in the nearest future. 
We also intend to use it to determine the distribution of 
the matter density (including dark matter) at the scales 
75/1-1 ^gj^g POTENT method. This hst was also 
used to estimate the cosmological parameters flm and 
CTg. The obtained constraint is given in graphical form 
on Fig. [71 The best numerical constraint is given by a 



combination Ss = {Qm/O-Sr^^as = 0.91 ± 0.05. 
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Fig. 10 Dependence of (^3-)"o"8 on a. The solid line de- 
picts constraint set by both models, the dashed line depicts 
constraint set only by the complex model only 
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Fig. 11 la confidence areas of S's vs. Qm dependence. The 
solid line corresponds to the complex model, the dashed line 
corresponds to the simple model 
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